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Abstract 

The paper proposes a simple and efficient method to study the dielectric properties of 
composite materials modeled using very large random resistor-capacitor (RC) networks. The 
algorithm is based on the Frank-Lobb reduction scheme and analyses the frequency dependent 
AC conductivity, permittivity and phase angle of composite materials. The simulation study is 
based on 1056 samples of random network containing 2 x 10 6 components randomly 
positioned in different proportion of capacitors. It has been found that these properties exhibit 
similarities to the universal dielectric response of dielectric materials. The results show that 
there is no variability between samples at low and high frequencies and across the power law. 


Key words: permittivity, network RC, conductivity, power law. 


1. Introduction 

Many materials used in daily life are composite and are often a blend of at least two 
constituents or phases [1] with conductive and dielectric (insulating) properties. These 
systems display a similar frequency dependent conductivity and permittivity [2, 3]. 

The frequency dependent complex conductivity a and permittivity e of composite materials 
are given respectively as [1, 8]: 


(5(m) = er(0) +A(im) n 

(1) 

£(co) = Sm n_1 + £ m 

(2) 


Where er(0) is the dc conductivity; A and B are constant; 0 < n > 1.0; £ m is the dielectric 
permittivity which exhibits a frequency independence at high frequencies, co is an angular 
frequency. 
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These electrical properties obey the anomalous power law [4] well known as the 
Universal Dielectric Response (UDR) [5, 6, 7]. Examples of universality of behavior 
of ceramic materials can be found in [2, 3]. Such features are modeled by simulating 
the electrical properties using a random network of resistors and capacitors [8, 9, 10]. 
The electrical network represents a heterogeneous microstructure that contains two 
constituents insulating and conductive regions [2]. An example of this mixture that 
can be shown as a resistor-capacitor network in Fig. 1(a). 

There are a large number of papers in the literature on the application of random RC 
network in modeling; however in most of those studies the authors use a small 
number of components (N) of networks and realizations. For example Almond et al. 
[14] explored a size of around 10 4 components, 100 realizations. Bouamrane et al. [4] 
used a network of 32 x 10 4 components, 256 realizations. Tuncer et al. [11] studied the 
response from several circuit cases with 26 xlO 4 lattice size and Hamou et al. [12] 
generated a network of 5 xlO 5 components, 10,000 realizations. An example of these 
results is shown in Fig. 2. The random network simulated consists of 512 components 
with 60%, 40% proportions of resistances and capacitances respectively. 


In this work, the authors have developed an algorithm to explore a large random RC 
network which includes 2 xlO 6 components and 1056 realizations. It has been found 
that for this size all the electrical properties of mixture materials are verified. The 
main goal of this work is to examine how well the response of conductivity, 
permittivity and phase angle of random mixtures materials may be modeled by large 
random network RC. 


2. The RC Network Model and the Calculation Process 
2.1 Principle of the numerical calculation 

A square network with N — 2 x n x n components is used to demonstrate the 
calculation methodology, where n is the number of nodes. For example, Fig. 1(b) 
shows a network with n = 5 and thus it has 50 components randomly positioned. The 
dielectric response £(u>) of a mixture random network R-C is determined by 
calculating the equivalent impedance or (admittance) from the following relation: 

£(u>) = £'(u>) + ie"(a)) — 1 /i£ 0 a)Z eq (co) — Y eq (co) /ie 0 co (3) 


The equivalent complex impedance Z eq (m) is calculated using an efficient algorithm, 
based on the Frank- Fobb (FF) reduction scheme. The FF method uses a mathematical 
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technique that simplifies the study of some electrical networks. This method consists 
of applying successive triangle-star Fig. 3(a-b) and star-triangle Fig. 3(c-d) 
transformations. 

The FL reduction scheme starts with the triangle of Fig. 3(a) which is transformed 
into a star to produce the configuration of Fig. 3(b) following the relation: (A — > Y) 


% 

^7 


* 1 * 2 +> 1 * 3 +> 2*3 
*3 

Y 1 V 2 +YiY 3 +Y 2 Y 3 

y 2 

y 1 y 2 +y x y 3 +y 2 y 3 
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( 4 ) 


Then the star of Fig. 3(c) is transformed back into a triangle as shown in Fig. 3(d) 
following the relation: (Y -> A) 


Y 9 = 
*io ' 

^ii = 


y 4 y 8 


y 4 +y s +y s 

Y s Y 8 

y 4 +y 5 +y 8 

y 4 y s 


Y.+Ya+Yf 
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( 5 ) 


The final result of this sequence of transformations is to reduce any large random 
network (RC) to a single link that have the same equivalent impedance as the entire 
network. 


2.2 Random network RC characteristics 

The square network (R-C) of Fig. 1(b) contains a random mixture of resistors and 
capacitors with a proportion of (1 - p) resistors ( 1 kD. ) and p capacitors (1 nF), with 

the admittance randomly as either y R -\/R or y c =icoC respectively where &>the 
angular frequency is 2nf . The critical percolation for a 2D square network is 
p c = 1/2 [14] and the fundamental study of random RC networks [3, 8, 9, 14] were: 

• When the proportion of capacitors is different from the critical proportion p =£ 
Pc'- 

i. If p < p c , the overall admittance is independent of m (|Y| <x y K ) at low 
and high frequency, the network is resistive. 

ii. If p > p c , the overall admittance is dependent on o) (|Y| oc y c ) at low 
and high frequency. 

• When the proportion of capacitors equal critical proportion (p = p c ): 
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i. At low frequency, half of the realizations of the network give an 
admittance response independent of co with a resistive percolation, and half 
gives an admittance response proportional to co with a capacitive 
percolation. A similar behavior is observed at high frequencies. 

• For all cases above, at intermediate values of co the admittances of capacitors 

and resistors become comparable( y H ~y c ), the emergent power-law behavior 

is observed and this is characterized by the admittance response |Y| oc co p . 

• For the resulting network associated with a simple logarithmic mixing rule, the 
complex conductivity and permittivity of the network could be described as: 

a* = {iooCy {R~ x f- p ex (im) p (6) 

s* = a*/ico oc (icu) p_1 oc (icu) p_1 (7) 

• The logarithmic mixing rule also predicts that the phase angle in radians, is 

00-0 = (1 -p).n/2 (8) 

Where p is the fraction of capacitors. 


3. Result and discussion 

To investigate the results obtained in previous work that uses small network sizes and 
realizations, the network responses of 1056 separate networks of randomly positioned 
resistors and capacitors were developed, which contains N =2097152 components. 

Fig. 4 shows the real £'(n>) and imaginary £"(n>) of permittivity as a function of the 
angular frequency co, expressed in log-scale, forp — p c — 1/2. Similar results are 
shown in Fig. 5 for p = 0.4 < p c ,p = 0.6 > p c . The logarithmic mixing rule of Eq. 
(6-7) used here leads to the power law frequency dependence on the complex 
conductivity and permittivity [8, 9]. These figures show an excellent dependence with 
co p and cjo~'° , where p and $ is the capacitor and resistor concentration respectively, 
as explained by Almond and al [9]. These results agree with those of the previous 
work. 

In Fig. 6(a), where the proportion is R =60 %, C — 40%, the phase of the 
conductive network at low and high frequencies is 0(m) * 0°, corresponding to the 
high and low frequency of admittance Fig. 7(a). At intermediate frequency 
(10 4 ~10 8 ) the power law of both s' and e" oc co -0,6 is preserved and Eq. (8) gives a 
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phase of —36°. The network simulated in Fig. 5 (b) exhibits no peak in dielectric loss 
s" because there is a conductivity of the conductive network. 

In Fig. 6(b), where the proportion is R — 40%, C — 60%, the phase is 0(n>) ~ —90° 
at high and low frequency, corresponding to the high and low frequency of s' in Fig. 
5(c) where the network becomes purely capacitive. At intermediate frequencies, the 
phase increases and is approximately constant in the frequency range (10 4 * * * ~10 8 ). In 
the same frequency range s' and e" lowering with frequency in power law s' and 
£"are proportional to m -0,4 . In this case Eq. (8) gives a phase 0(m) = —54°. For the 
case of R=50%, C=50%, as shown in Fig. 6(c), at low frequency the phase of the 
network is approximately —90° therefore it is influenced by the high impedances of 
the capacitors. At the high frequencies the phase approximately 0°, and the network 
is influenced by the high impedance of the resistors. At intermediate frequencies the 
phase equal ~45° Fig. 6(c), predicted by the Eq. (8) and both of s' and s" go with 
power law proportional to m -0,5 Fig. 4. 

The results of the admittance in cases where the proportions of capacitors p — 0.4, 
p = 0.6, p = 0.5 are shown in Fig.7 for 1056 different realizations of large network 
size of N — 2097152 components as a function ofrn. By plotting the numerical 
calculations with the logarithmic mixing rule of Eq. (6), it can be seen that the 
predictions of this rule perfectly agree with results of the numerical simulations. In all 
cases, the intermediate frequency power law region is clearly defined with o) p . 


4. Conclusions 

This work has analyzed very large number of random network resistors -capacitors 
and the frequency responses of the permittivity, admittance, and phase angle are 
simulated. The calculations have been performed using an efficient algorithm 

developed based on the Frank-Lobb reduction scheme. The results of this large 

network confirm that the power law characteristic is an emergent property of random 

RC network and the frequency response of permittivity, admittance, and phase angle 

of composite materials are verified. It has been shown that the permittivity decreases 
with m - ’ 9 , d is the proportion of resistance and the admittance increase with cj p 
where p is the proportion id capacitance. It has been found that frequency ranges 
through which these power law regions are produced do not increase with the size of 
the network. The prediction of the logarithmic mixing rule is clearly demonstrated for 
large random networks RC. For these reasons it can be concluded that the proposed 
method is valid for networks of any sizes. 


5 



Acknowledgments 

M.Aouaichia and R. Bouamrane would like to thank the Directorate General for 
Scientific Research and Technology Development (DG-RSDT) for providing our 
university with a high performance computing facility. 


References 

[1] E. Tuncer, Y.V. Serdyuk, S.M. Gubanski, Dielectric Mixtures: Electrical 
properties and modeling, IEE Transactions on Dielectrics and Electrical 
Insulation. Vol. 9 No. 5 (2002). 

[2] C. R. Bowen and D. P. Almond, Modelling the ‘universal’ dielectric response in 
heterogeneous materials using microstructural electrical networks. Materials Science 
and Technology. Vol. 22 No. 6 (2006). 

[3] S. PANTENY, R. STEVENS, AND C. R. BOWEN, the Frequency Dependent 
Permittivity and AC Conductivity of Random Electrical Networks. 319 (2003) 
199-208. 

[4] R. Bouamrane and D. P. Almond, The ‘emergent scaling’ phenomenon and the 
dielectric properties of random resistor-capacitor networks J. Phys. Condens. Matter, 
(2003), 15, 4089-4100. 

[5] A. K. Jonscher, The “Universal” Dielectric Response: Part I. IEEE Electrical 
Insulation Magazine March/April (1990) Vol. 6 No. 2. 

[6] K.L.Ngai, A.K.Jonscher, C.T.White, On the origin of the universal dielectric 
response in condensed matter. Nature Vol. 277 (1979). 

[7] A.K. Jonscher, Universal Relaxation Law: A Sequel to Dielectric Relaxation in 
Solids, Chelsea Dielectrics Press, (1996). 

[8] B. Vainas, D.P. Almond, J. Luo, R. Stevens, An evaluation of random R-C 
networks for modelling the bulk ac electrical response of ionic conductors, Solid State 
Ion. 126 (1) (1999) 65-80. 

[9] D.P. Almond, B. Vainas, The dielectric properties of random RC networks as an 
explanation of the ‘universal’ power law dielectric response of solids, J. Phys.: 
Condens. Matter 11 (1999) 9081-9093. 

[10] Jeppe C. Dyre and Thomas B. Schrpder, Universality of ac conduction in 
disordered solids, Reviews of Modern Physics, Vol. 72 No. 3 (2000). 

[11] Enis Tuncer, dielectric mixtures importance and theoretical approaches, IEEE 
Electrical Insulation Magazine Vol. 29 Issue 6. 


6 



[12] R.F. Hamou, J.R. Macdonald, and E. Tuncer, \Dispersive dielectric and 
conductive effects in 2dr esistor-capacitor networks," Journal of Physics: Condensed 
Matter, vol. 21, no. 2, p. 025904 (13pp), 2009. 

[13] D.J. Frank and C.J. Lobb. Highly efficient algorithm for percolative transport 
studies in two dimensions. Physical Review B , 37(1) (1988) 302-307. 

[14] D.P. Almond, C.J. Budd, M.A. Freitag, G.W. Hunt, N.J. McCullen, N.D. Smith. 
The origin of power-law emergent scaling in large binary networks. Physica A 392 
(2013) 1004-1027 


7 



Figure 



(b) 




X*T 


n- 1 I 


O it 


7 II 


X 


II — — It 


X 


I — • 


Fig. 1. (a) Microstructure of Al 2 0 2 —Ti0 2 composite: grey phase corresponds to ( C AL202 , capacitor); 

white phase corresponds to ( Ti0 2 , conductor) [2]. (b) Example of resistor (white phase) capacitors (grey 
phase) network. 
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Fig. 2. Example of Ac conductivities of 256 R-C networks of 512 components (60% R, 40% 
C). 
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Fig.3. Fundamentals transformation: (a — b) -* (A — Y), (c — d) -> (Y — A). 



Figure 


p = 0.5 




Fig.4. Real and imaginary permittivity £'(n>), £"(n>) as a function of m with rresponses of network 
simulation of p = 1/2. The results of Fig. 4 are obtained from the simulation of a network 
containing 2097152 components for 1056 random realizations. The green line shows the 
logarithmic mixing formula, Eqn (7), prediction for the networks. 
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Fig. 5. The real and imaginary permittivity £'(cu), £"(m) as a function of (D with 
Responses of network simulation of p — 0.6, 0.4. The results of Fig. 5 are obtained 
from the simulation of a network containing 2097152 components for 1056 random 
realizations. The green line shows the logarithmic mixing formula, Equation (7), 
prediction for the networks. 
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Fig. 6. The phase-frequency of network which contains 2097152 components, with response of 
1056 realizations for different cases of p — 0.4, 0.6, 0.5 respectively. The green line shows the 
logarithmic mixing formula, Eqn (8), prediction for the networks. 
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Fig. 7. The overall admittance of network as a function of frequency which contain 
2097152 components, with response of 1056 realizations for response of p — 
0.4, 0.6, 0.5. The green line shows the logarithmic mixing formula, Equation (6), 
prediction for the networks. 


